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Abstract 

We present the full numerical solution for the 15- dimensional space of linearized de- 
formations of the Klebanov-Strassler background which preserve the SU (2) x SU(2) x Z2 
symmetries. We identify within this space the solution corresponding to anti-D3 branes, 
(modulo the presence of a certain "subleading" singularity in the infrared). All the 15 in- 
tegration constants of this solution are fixed in terms of the number of anti-D3 branes, and 
the solution differs in the UV from the supersymmetric solution into which it is supposed 
to decay by a mode corresponding to a rescaling of the field theory coordinates. Deciding 
whether two solutions that differ in the UV by a rescaling mode are dual to the same theory 
is involved even for supersymmetric Klebanov-Strassler solutions, and we explain in detail 
some of the subtleties associated to this. 
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1 Introduction 



Antibranes in warped deformed conifold Klebanov-Strassler (KS) backgrounds [1] are a staple 
ingredient of string phenomenology and cosmology constructions, being essentially the only 
method for lifting AdS solutions with stabilized moduli, to dS solutions, and thus give rise to a 
landscape of dS vacua of string theory [2]. 

Over the past few years we have undertaken a programme to construct the full space of 
first-order SU(2) x SU(2) x Z2-invariant deformations around the KS background, in order to 
establish whether a solution corresponding to anti-D3 branes in this background exists, whether 
it has the properties one expects from the brane-probe analysis of [3], and whether it is dual to a 
metastable vacuum of the dual boundary theory The underlying philosophy of this programme 
has been that one cannot decide a-priori that a metastable anti-D3 brane solution must exist, 
and then accept whatever boundary conditions are necessary in order for this to happen, but 
rather one should start from a set of physical infrared and ultraviolet boundary conditions, and 
ask whether a solution compatible with these boundary conditions exists or not. 

The key results of this investigation have been: 

1. One can find all the homogeneous solutions to, and thus solve implicitly the equations [4] 
governing the first-order perturbations. The full solution seems at first to involve 8 nested 
integrals [5]. 

2. One can simplify these and write the full solution in terms of 2 nested integrals [6], which are 
in fact integrals of rational functions multiplying the warp factor and Green's function of the 
KS background. 

3. One can write the UV and IR expansions of the generic solution to this space of deformations, 
and identify all the UV normalizable and non-normalizable modes, as well as the infrared physical 
boundary conditions for D-branes [5]. 

4. The force on a probe D3 brane in the first-order perturbed background depends only on one of 
the 16 integration constants, and this constant must be nonzero if the solution is to correspond 
to antibranes [5]. Furthermore, the full functional expression of this force can be calculated [7], 
and matches exactly the expression one obtains from "Newton's third law" arguments a la 
KKLMMT [8]. 

5. The putative solution for anti-D3 branes smeared on the three-sphere at the tip of the 
KS solution is expected to have a singularity in the five-form and warp factor, coming from 
the physical brane sources. Besides this, the solution must also have a subleading singularity, 
proportional to the coefficient of the brane-attracting mode of the solution. 

As explained in [5], if the singularity is not physical, then the backreaction of anti-D3 branes 
in the KS solution gives rise to a large deformation of this solution, which cannot be captured in 
perturbation theory, much like when one tries to construct metastable vacua using type IIA brane 
engineering [9]. On the other hand, if the singularity is physical, then our technology produces 
the full first-order backreacted solution corresponding to antibranes in the KS background, as 
well as all first-order deformation of the KS solution by non-normalizable SU (2) x SU (2) x Z2- 
invariant modes, corresponding to all the relevant and irrelevant deformations of the dual field 
theory. 

This subleading singularity cannot be attributed to any brane source (it has the wrong 
orientation), or to brane-flux annihilation (it is linear in the antibrane number, while the brane- 
flux annihilation is nonlinear). However, as mentioned in [5] and argued in [10], it is possible 



2 



that this singularity is an artifact of perturbation theory, and may not be present in a fully- 
backreacted solution for antibranes. On the other hand, obtaining a fully-backreacted solution 
for antibranes in ISD flux backgrounds seems to run into trouble in less complicated setups 
[11, 12], and can even be ruled out by topological arguments (that yield a physics similar to 
the one found in [9]). If the results of [12] extend to the KS solution, then the presence of a 
subleading singularity in perturbation theory will look with hindsight as an indication of a more 
profound problem with the whole construction. 

Given that the arguments about this singularity fall mainly outside of the scope of our 
perturbation theory machinery, it is best to hedge our bets both ways, and ask whether inside 
the 15-dimensional space of parameters that characterize our first-order solution one can identify 
a solution that has the correct physics to correspond to anti-D3 branes in the KS geometry, 
subleading singularity aside. Identifying this solution inside the 15-dimensional space is simpler 
than finding a needle in a haystack, but not by far: One has to throw away divergent terms both 
in the UV and in the IR expansion [5], and to impose the correct D-brane boundary conditions 
on the divergence of the warp factor and electric field in the infrared. 

Those conditions yield algebraic relations between the various integration constants that 
appear in the UV or IR expansions of the fields; however, the integration constant that appears 
in the UV expansion of a given field, say the dilaton, is not the same as the one that appears in 
its IR expansion, but differs by highly nontrivial combination of the other integration constants. 
Hence, even if we impose all the physical boundary conditions in the UV and in the IR, we 
are far from being done, because the UV conditions are expressed using the UV integration 
constants, and the IR conditions are expressed using the IR integration constants, and it is 
possible that upon translating the UV conditions into IR variables one may have the unpleasant 
surprise that these conditions are incompatible. Hence, in order to identify the correct antibrane 
solution inside the 15-dimensional space of first-order deformations, and to establish whether 
this solution is dual or not to a metastable vacuum of a supersymmetric field theory, it is crucial 
to relate the UV and IR expansion coefficients, which is the main purpose of this paper. 

Before unveiling those results, we would like to point out that identifying whether two 
asymptotically-KS supergravity solutions are dual to vacua of the same field theory is not as 
straightforward as it might seem, even for supersymmetric solutions, essentially because, besides 
the seven normalizable and seven non-normalizable deformations, there exists another deforma- 
tion corresponding to rescaling the field theory directions. Of course, if two solutions differ by 
non-normalizable deformations, they clearly are dual to two different field theories; however, as 
we will explain in Section 5, two solutions with different rescaling parameters may or may not 
belong to the same theory. Hence, using purely UV data one cannot distinguish asymptotically- 
KS supersymmetric solutions that we expect [13] to be dual to different field theories, unless one 
introduces extra assumptions about the infrared of the solutions, or about their bulk behavior. 

Anticipating our results, we compute the unique solution that has the correct infrared and 
ultraviolet divergences (modulo the subleading singularity) to describe anti-D3 branes in the KS 
background. All the parameters of this solutions can be determined in terms of the number of 
antibranes. Nevertheless, much like for supersymmetric solutions, one cannot distinguish using 
purely UV data whether this solution describes a non-supersymmetric vacuum of a supersym- 
metric solution, or whether it is dual to a distinct non-supersymmetric theory. To achieve that 
one must therefore introduce extra assumptions about the infrared or about the bulk. 

In Section 2 we give a lightning review of the general construction of first-order deformations 
around the KS solution (the full details can be found in Appendix A), and in Section 3 we review 
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the simplified analytic solution found in [6] in terms of two nested integrals (whose full details can 
be found in Appendix B). In Section 4 we explain the procedure we use to relate the UV and the 
IR integration constants, and illustrate with more details how this procedure can be implemented 
for one of the perturbation modes. We also give the relations between the UV and IR integration 
constants of the other modes; the derivation of all these relations is left for Appendix C. In 
Section 5 we present the different criteria for distinguishing supersymmetric asymptotically- 
KS solutions, and in Section 6 we identify the solution for anti-D3 branes inside the space of 
solutions. Section 7 is devoted to the relation between our solution and the one obtained in [14] 
by perturbing around the Klebanov-Tseytlin (KT) solution, and to the identification within our 
space of solutions to perturbation of the KS solution by non-normalizable modes dual to gaugino 
masses. 



2 Non-supersymmetric deformations around the 
Klebanov-Strassler background 

2.1 Ansatz and background solution 

We wish to construct the backreacted solution corresponding to N anti-D3 branes smeared on 
the S 3 at the tip of the warped deformed conifold. We use the Ansatz proposed by Papadopoulos 
and Tseytlin [15], which is the most general one (with vanishing RR axion Co) that preserves 
the SU(2) x SU(2) x Z2-symmetry of the Klebanov-Strassler solution (KS). The metric is 

ds\ Q = e 2 ^"* dsl 3 + e" 6 ^ rfr 2 + e* + * {g{ + gj) + e— (<£ + gj) + g\ , (1) 

where all the functions depend on the radial variable r. The fluxes and the dilaton are 

#3 = \ (k - f) g 5 A ( gi Ag 3 + g 2 A gt) + dr A (f gi Ag 2 + k'g 3 A g 4 ) , 
F 3 = Fg 1 Ag 2 Ag 5 + (2P-F)g 3 Ag 4 Ag 5 + F'dTA(g 1 Ag 3 + g 2 Ag i ), (2) 

"vrQ 



F 5 = J- 5 + *J5, ^5 

$ = $(r), C = 0, 
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+ (k-f)F + 2Pf 



9iAg 2 Ag 3 Ag 4 Ag 5 , 



where P, Q are constants while /, k and F are functions of r. A prime denotes a derivative with 
respect to r. 

The fields from this Ansatz are collectively denoted a , a = 1, 8. We will study and fully 
determine the solution space of first-order non-supersymmetric deformations of the supersym- 
metric Klebanov-Strassler theory, 

r = r + r^z) + o(z 2 ) . (3) 

The background fields 0q are given by the Klebanov-Strassler solution without mobile D3- 
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branes: 



e x ° = \Kt) 1/2 (± sinh(2r) -r) V3 
e yo = tanh(r/2), 

sinh(2r) -r) V3 

e 6po _ 24- 



e 



h(r) sinh r 

6Ao " ^Mt) (|sinh(2 r)-r) 2/3 sinh 2 r, (4) 
(r cothr — 1) (coshr — 1) 



h = -p 

k = ~P 



sinhr 

(r cothr — 1) (coshr + 1) 



F = P 



sinh r 
(sinh r — r) 

sinh r 



$o = 0, 

g = o, 

where Eq is the deformation parameter of the conifold, related to the confinement scale of the 
dual gauge theory. Of significance are also the warp factor h and the Green's function j for this 
background: 

h{r) = 32 P 2 f g (cosh u sinh u — u) 1 ^ 3 du , (5) 

J T sinh u 

Note that the last equality in (4) implies we are taking g s — 1. Furthermore, the dimensionful 
constant P is related to the quantized dimensionless units of flux M entering in the rank of the 
gauge groups of the dual field theory (see section 5.2) by 

P = \Mol, (7) 

So as to avoid extra clutter, in what follows we take a' = 1, and £q = 1. 



2.2 First-order perturbation equations, conditions and physical sig- 
nificance of the integration constants 

Using a method due to Borokhov and Gubser [4] and reviewed in the Appendix A, finding 
linearized deformations away from a supersymmetric solution, can be reduced to solving two 
sets of first-order ordinary differential equations in the radial variable r, instead of second-order 
differential equations. Out of those two sets, the first one forms a closed system for the variables 
£ a that can be thought of as "conjugate momenta" for the perturbations 0" of the fields entering 
our Ansatz (1), (2). The integration constants associated to that first system are labelled X a , 
and are non-zero for a non-supersymmetric solution. The integrations constants from the second 
system of coupled lst-order ODE's are denoted Y a . 
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For the problem of present interest, i.e. the backreaction of anti-D3's on KS, the solution to 
the system was found in [5, 6] after applying the following change of basis 1 

x - 2p - 5 A, y, x + 3p, x - 2p - 2 A, f + k + F, $ 

There is one relation between the constants X a that has to be obeyed on the whole space of 
solutions. Namely, the zero-energy condition 

6X 2 -4X 3 -6PX 5 -9PX 7 = . (9) 

Another integration constant, Y\ as it happens, looks naively like it can be gauged away by a 
rescaling of the four- dimensional coordinates but as we will see later plays a crucial role in the 
physics. We are therefore left with fifteen meaningful integration constants. 

Out of those fifteen parameters, the one called X\ plays a key role. Indeed, the force exterted 
on a probe D3-brane is directly proportional to it and does not depend on any other integration 
constant [5]. Its expression was found in [7] and is given by 

2 e~ 2x ° 

I e~ 2 ^ X 1 h{r) , 

32 ^ (10) 

3 (sinh2r-2r) 2 / 3 ' 1 ; 



'D3 + = 36-^1 



One can also use the conventions of [16] to describe the same result for a first order expansion 
around any warped Calabi-Yau background with ISD flux. Here the derivative of the DBI and 
WZ actions for D3-branes are respectively proportional to the warp factor e 4A and the four- form 
RR potential C4 = adx° A ... A dx 3 , where in the language of (1) and (2), we have 



A = A + p--, a ' = -e 4A+4p ~ 4x 
2 



^R + kF + f (2P-F) 



The force is found to be 

F m± = $' T , where $± = e 4A ± a , (12) 

and by D3_ we mean D3-branes. The combinations $± are sourced by D3± respectively, and 
by |G ± | 2 [17, 10] where G± = G 3 T % * G 3 and G 3 = F 3 + ie^H^ 



3 Our analytic solution for the full space of first-order 
deformations around KS 

In a previous work [6] , we found that the fully analytic generic solution to the most general first- 
order deformation of the Klebanov-Strassler background involves at most two nested integrals 
of the form 

/T PT 
h(u)f{u)du, or / j(u)f(u)du, (13) 
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Figure 1: The profile of the field 8 corresponding to a shift of the dilaton, for the following 
choices of integration constants (with e.g. P = 1). Blue, also labelled (a): X\ = 1,X 5 = 
-f ,X 6 = X 7 = 5,X 8 = 2, Y 8 = -88.05; Red ®: X 1 = X 6 = X 7 = 1,X 5 = -f,X 8 = 1.8, Y 8 = 
-111.5; Yellow ©: Xi = X 7 = 2,X 5 = -|,X 6 = 8.608, X 8 = -0.843, Y 8 = -133.9. In each 
case, Y 8 is fixed so as to ensure that 8 (oo) = 0. 



where f(r) is a certain combination of hyperbolic functions. Expressions for the warp factor 
h(r) of the KS background and its Green's function j(r) are provided in (5) and (6). 

Let us illustrate this with the result for 8 , corresponding to shifts in the dilaton. The 
analytic solutions for all seven remaining modes are consigned to Appendix B, and more details 
of the derivation can be found in [5, 6]. 

4> 8 =Y 8 -64X 8 j( T ) + ^-h(r) 

„ ,, f T (u cothu — 1) 

-64PX 6 / — ^ — du 

J i sinh u (coshw sinhw — u) 1 

+ 2 — j h(r) H Xi csch 2 r (cosh r sinh r — r) 1 ^ 3 h(r) 

P 3 

64 , s , s 32 /" T (sinh 2 M + 1 - u cothw) 

+ -X 1 / i (r) J (r)--X 1 / / — ^s/i(«)d«. (W) 

y J Ji sinh m (cosh-u sinnu- u) 

We have chosen to integrate in the domain [l,r], given that many of the integrands (like the 
one from the last term above) are infrared- divergent. Once the limits of integration are fixed, 
the constant Y 8 in (14) is defined unambiguously. The profile for 8 is given in Figure 1. 

The infrared and ultraviolet behaviors of the modes are given in Appendix C. Some of the 
integration constants appearing in the infrared expansions (like Y^ R or Y 6 /i? ) correspond to 



1 Note that (j>4 is the perturbation to the warp factor, namely (f>4 — —2 A, since the warp factor of the KS 
theory (5) is such that h(r) = e 4A — e 4 -4+4p-2a;^ Q£ a j g0 e q ua tion (11) below. 
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unphysical divergences of various fields, and we will set them to zero. Other constants (like 
Y 7 IR or Xi) correspond to physical divergences in the warp factor and in the RR five-form field 
strength coming from the presence of smeared anti-D3 branes, and we need to keep them in 
the final solution. We will explain this procedure when we construct the antibrane solution in 
section 6. 

In order to stress out how the integration constants X a and Y a are paired into normalizable 
and non-normalizable modes we also remind the reader of the UV behaviors of those modes [5], 
which one can also extract from the expansions in Appendix C: 



dim A 


non-norm / norm 


integration constants 


8 


ry^ j ry 8 


Y A /X l 


7 




Y 5 /X 6 


6 




X 3 /Y 3 


5 


ry 1 ' ry 5 




4 


y»0 / ry 4 


Y7, Yg, Y 1 / A 5 , X4, X 8 


3 


ry 1 / ry 3 


X2, X 7 /Yq, y 2 


2 


ry 2 / \y- 2 





Table 1: The UV behavior of all sixteen modes for the SU(2) x SU(2) x ^-symmetric deformation 
Ansatz around the Klebanov-Strassler solution. 



4 Relating the IR and UV integration constants 

Given that ultimately we will have to impose boundary conditions on the generic analytic so- 
lution to the full space of first order deformations around KS, we should look at the IR and 
UV behavior of the modes <fi a . Their somewhat lengthy analytic expressions are gathered in 
Appendix B and were first found in [6]. Moreover, it is not enough to consider the expansions 
shown in Appendix C. The zeroth-order terms in the expansions collected in that Appendix in- 
clude arbitrary integration constants coming from indefinite integrations, which are generically 
denoted as Y^ R , Y^ v ■ In order to determine how the Yj R 's are related to the Y^ v, s and thus to 
connect the IR and UV regions, we have to perform a numerical integration that will fix Y^ v 
as follows: 

Y a UV = Y' R + iTN a b X b , (15) 

6=1 

where N is a matrix of numerical coefficients arising out of evalutions of the single and double 
integrals appearing in the analytic solutions for the a modes. 

4.1 Our results 

All in all, following the procedure we have just outlined, the relations between all 2 the Y^ v and 
Y* R that we have derived are as follows: 

2 Except Y4, which is far more difficult to get and will not be needed for our following analysis in any case. 
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{y 8 uv \ 

vuv 
vuv 

yUV 

xrUV 
J 5 

yUV 



XrlR 
*2 



\ 



VlR 



2Yj* 

5 V IR 
3 r 2 



ylR 
. r 7 



yIR + 

+ " " 



— P Yj R + P YJ R 



3_P V IR 
2 J 2 



I 8 



+ N- 



x 2 
x 3 
x 4 
x 5 
x Q 
x 7 

W 



(16) 



with the matrix N 
N = 



/ -235.3 P 2 
-3.870 P 2 
93.63 P 2 
-123.8 P 2 
-165.9 P 3 
100.6 P 3 
y-225.8P 3 






250.0 
-40.25 
-20.16 P 
-166. 7P 
83.34 P 









-36.47 P 


35.71 P 


-18.24P 


83.34 


7.791 


83.34 P 


-12.37P 


166. 7P 


206.7 


93.84 


-284.0P 


61.22P 


-243.8 P 


70.31 


-1.827 


22.93 P 


35.50 P 


71.33 P 


19.52P 


1.488 P 


14.08 P 2 


11.90 P 2 


36.32 P 2 


81.27P 


-46.06 P 


221.4 P 2 


-48.57 P 2 


265.8 P 2 


-94.65 P 


16.52 P 


-158.9 P 2 


35.92 P 2 


-221.4 P 2 



53.56 \ 




17.85 P 
-8.545 P 
17.09 P 



/ 



The above relations (16) depend at an intermediary stage on our results for the relation 
between the integration constants Y a that appear in the analytic solution (131)-(143) and the 
constants Y^ R that appear in the IR expansions (146)-(153), obtained via the method summa- 
rized at the beginning of this section and further expanded upon in the next subsection. We 
provide them here as a matter of having accessible intermediate results: 

X 2 
X 3 

x, 
x 5 
x 6 
x 7 

w 



/y 8 ir \ 


( Ys \ 


vIR 
I 2 


Y 2 


XrlR 
J 3 


2Y 2 + Y 3 


ylR = 


Y x 




Y 5 - f Y 8 


\rIR 


-%Y 2 + Y 6 


Vr IR ) 


\PY 2 + Y 7 ) 



+ M 



(Y'R.Y) 



(17) 



9 





352.6 P 2 













25.86 P 2 





-33.23 


3.918 




-18.62 P 2 


-99.69 


15.54 


-0.9673 




144.4 P 2 


98.79 


-67.47 


5.146 




92.62 P 3 


12.26 P 


-9.501 P 


-4.435 P 




8.129P 3 


24.44 P 


-1.632 P 


1.147 P 




-1.307 P 3 


-38.81 P 


-4.754 P 


3.491 P 



Analogously, the link between the parameters Y a 
UV/IR relation (16). 



36.47P 


-41.56 P 


18.24 P 


-53.56 


-38.81 P 


-3.432 P 


-69.25 P 





7.797 P 


-7.959 P 


15.92P 





-81.34 P 


-44.35 P 


-153.9 P 





-16.54 P 2 


-18.03 P 2 


-22.52 P 2 


-11.85 P 


-4.773 P 2 


2.180 P 2 


-11.20 P 2 


-3.979 P 


1.749 P 2 


3.599 P 2 


-6.256 P 2 


7.959 P 



and Y a can similarly be obtained from the 



4.2 An illustration of the procedure 

As an example making this procedure plainer to the reader, we show how we relate Y% v and 
Yg /R . This is a three-stage procedure: 

(i) first, we relate Y% R and the parameter Y 8 appearing in (14); 

(ii) we next obtain the relation between Y% and Y$; 

(iii) finally, using results from the above steps, we get 3g in terms of Y$ . 

In order to implement step (i) above and relate Y$ R to Y 8 , we expand the integrands entering 
the IR expansion of the solution to the 8 equation up to a certain power in r. We then evaluate 
the indefinite integral and call Y$ R the constant term in <j>&. The first few terms in those 
expansions are given by (146), which we provide here for convenience: 

% R = \ (f (f) (3PX 6 - hoX,) + 32 ■ 2 1 / 3 . 3 2 / 3 X 8 ) + Yj* + 0{r) . (18) 

We now have to match (18) at some small r with the numerical value of (f) 8 that we obtain by 
performing the integrals in (14) numerically. Since the expansions for the integrands are good 
up to t > 1, we did choose to match at r = 1, where the integrals that enter the solutions for 
the 0's are zero by definition. Evaluating numerically (14) at r = 1, we find 

8 (r = 1) = Y 8 + 84.0493 P 2 X l + 28.5159 PX 5 + 14.2579 P X 7 + 41.2221 X 8 , (19) 
while from the IR expansion of <p 8 (18), we have 

0^(1) = Y^ R - 268.524 P 2 X 1 - 7.9588 PX 5 + 41. 5621 PX 6 

-3.97940 PX 7 + 94.786 X 8 . (20) 

Comparing the above two results, (19) and (20), we finally obtain the end-result of step (i) 
above: 

Y 8 IR = Y 8 + 352.574 P 2 X 1 + 36.4747 PX 5 - 41.5621 PX 6 

+ 18.2373PX 7 - 53.5642X 8 . (21) 

With this relation at hand, we can furthermore make sure, as one more consistency test, that 
the numerical integrals and the series agree at small r. The result is shown on Figure 2. 
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We go through the same recipe for the UV and compare the value of the UV series of the 
integrands with the value of 0§ that we have obtained by performing the integrals numerically 3 
at r = 15. When the dust settles down, we find the following relation between Y% v and Y s : 

Y% v = Y 8 + 117.318 P 2 X l - 5.85263 P X 6 . (22) 

As one extra check, inserting the above result in the UV expansions, we can verify that the 
UV series approximates well our numerical results at large r. This can also be see on Figure 2. 

Note that for 0§ there is a rather large range of overlap between its IR and UV series 
expansions. So, with hindsight, for this particular mode, we could have avoided going through 
tedious numerical work. On the other hand, for most of the other a fields, the overlap is 
much narrower. Therefore, in order to attain satisfactory precision in relating the IR and UV 
integration constants, we have opted for a careful numerical analysis. 



08 (r) 




-100 L 

Figure 2: The numerical solution for the field (f>$ for X x = 1,X 5 = — X 6 = 5,X 7 = 5,X S = 
2, Y 8 = — 88.05, P = 1 (underlying blue solid line). The red and orange dashed lines correspond 
respectively to the IR and UV expansions. 



5 Asymptotically KS solutions and their field theory in- 
terpretation 

Having found the full 15-dimensional space of perturbative solutions around the KS background, 
we would now like to develop the machinery that will allow us to identify whether the antibrane 
solution is in the same theory as the supersymmetric background into which it is conjectured to 
decay [3]. However, as mentioned in the introduction, distinguishing between asymptotically- 
KS solutions and arguing which background is dual to which field theory using only UV data 
is not trivial even for supersymmetric solutions, essentially because of the existence of the scale 

3 With as much precision as desired. Here, for both IR and UV expansions, we have settled for 20 orders of 
WorkingPrecision using Mathematica. The UV series expansions were derived up to order 15. 
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deformation Y\_, which equivalently can be traded for the e parameter that characterizes the size 
of the deformed conifold before the warping. 

If two solutions differ by non-normalizable deformations, they are dual to two different field 
theories. However, our fifteen-dimensional deformation space has the peculiarity that there are 
seven pairs of normalizable/non-normalizable modes and then one extra mode Y±. The putative 
partner to Y\ is eliminated by the zero-energy condition and it may seem that Y\ itself is a gauge 
artifact which can be removed by rescaling the four-dimensional space-time coordinates. As we 
will mention in more detail below, while for a single vacuum this is true, if there are two isolated 
vacua in the same theory then there remains a dimensionless number (essentially the ratio of 
the confinement scales) which can be attributed to Y\. 

One can inquire whether two solutions that have the same non-normalizable modes but two 
different e's, hence two different scale deformations, are dual to the same field theory. The 
answer is not clear, because one can change e and at the same time change also the number of 
mobile branes, keeping the total charge at infinity constant. Changing e changes the volume of 
the space, and since the space has charge dissolved in flux, one also changes the total charge; 
one can compensate for this change by introducing or taking away mobile branes. 

Hence, a vacuum with no mobile branes and one value for e has exactly the same UV data 
as a vacuum with one mobile brane and another value of e, or a vacuum with, say, 17 mobile 
branes and yet another value of e. Clearly these solutions cannot be all dual to vacua of the 
same KS field theory. On the other hand, a background with M mobile branes (where M is 
the amount of RR three-form flux on the KS three-cycle) and a certain value of e and another 
one with no mobile branes were argued in [13] to be dual respectively to the mesonic vacuum 
and the baryonic vacuum of the same SU(kM) x SU (kM + M) theory. Hence, even in the 
supersymmetric theory, one cannot decide whether two vacua with different scales and different 
amounts of mobile branes are in the same theory by simply examining their UV data. 

In this section, we discuss the supersymmetric KS situation in detail, and argue that in 
order to be able to use UV data to distinguish between two supersymmetric asymptotically-KS 
solutions that should not be dual to the same theory, one must introduce an additional criterion. 
The most obvious choice is requiring that the value of the NSNS B 2 field that wraps the S 2 which 
shrinks to zero size at the conifold tip must be zero, and can only jump by integral periods. After 
all, the S 2 is topologically trivial, and if the integral of B 2 is nonzero, one can stay at a fixed 
radius, consider a very small closed fundamental string at the north pole and take it around the 
S 2 to the south pole; during this process its world-sheet action will pick up a phase proportional 
to the B 2 integral. If one now brings back the string to the north pole, the string will interfere 
destructively with itself unless the integral of B 2 on S 2 is an integer 4 . This argument is similar to 
that ruling out Dirac strings, and in principle should also hold in the presence of D3 or anti-D3 
branes. 

A second possible criterion is requiring that the integrals of the H 3 from the origin to a 
certain holographic screen differ by an integer amount for two solutions in the same theory, or 
equivalently that the difference in the number of Seiberg duality cascades between two solutions 
dual to vacua of the same theory has to be integer-valued. This criterion has a clear physical 
justification for compact settings, where the KS throat is seen as the zoom-in of a compact CY, 
and where the three-cycle wrapped by H 3 that appears non-compact from a KS perspective is 
in fact embedded into a compact CY three-cycle. However, for a non-compact KS solution this 
criterion is very hard to justify from a holographic perspective, because it involves integrals over 

4 We thank Nick Warner for this argument. 
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the whole bulk. 

One can also use the analysis of [13] to reverse-engineer a criterion that allows one to distin- 
guish between vacua with various numbers of mobile branes without introducing any extra IR 
boundary conditions, and using only UV data. This third criterion (summarized in equation (53) 
for the first mesonic vacuum), if correct, certainly requires a more physical explanation. 

Of course, another possibility is that the holography is just not refined-enough to distinguish 
between these different theories, especially because we are dealing with cascading solutions that 
are not asymptotically AdS, cannot be thought of as the near-horizon of any brane, and have 
an infinite charge unless one imposes an UV cut-off. 

In this section, we will use the first criterion, and give a holographic recipe for distinguishing 
between asymptotically-KS vacua that have different numbers of mobile branes. 



5.1 Maxwell charge, Page charge and mobile D3-branes 

For a supergravity solution with non-trivial Wess-Zumino terms one can generally define three 
different types of charges [18, 19], which we review in this section. The D3-Page charge, special- 
ized to the KS background is 



(23) 



This is conserved and is independent of the radius at which it is evaluated. In string theory it 
must also be quantized. If we shift B 2 by a small gauge transformation B 2 — > B 2 + dh\ for some 
one-form A 1; the charge stays invariant. In principle there are two independent ways to generate 
a non-zero, integer-valued starting from the smooth KS background: 

^5 -> ^5 + 27Qttvo1 t i,i , (24) 
B 2 -> B 2 + ^iru 2 , (25) 

=► QdT = Q-P (26) 

where (Q,p) G Z 2 , M is related to P by (7) and 

vol T M =^ gi A 02 A 5-3 A g 4 A g 5 

u 2 =\i.9i A #2 + #3 A #4) • (27) 

Having Q ^ generates a singularity in both the warp factor and *J~s, which one must 
interpret as due to Q D3 branes smeared on the tip of the deformed conifold. On the other 
hand, the meaning of the singularity due to p 7^ is more subtle, and if one imposes as an IR 
regularity condition that the B 2 field at the KS tip be zero or an integer mod M, then Q P D l ge = Q 
measures the number (modulo M) of mobile BPS D3-branes in any particular KS background. 

The Maxwell D3-charge is 

Q%S" = 7iLt (28) 



(4vr 



where the integral is performed on a Gaussian surface at the UV cut-off r = r c . There are two 
physically distinct contributions to the Maxwell charge, from mobile branes (g&) and from charge 
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dissolved in flux (qf): 



QMa* = qb + qfj (29 ) 
1 

(47 2 



ft = 77"^ / -- F ^ ( 3 °) 



* = (^(i/ 5 -4^ = (^l/ 3AF3 - (31) 

The Maxwell charge depends on the scale at which it is measured, but if we fix a holographic 
screen, we expect physical processes to preserve its value at the screen. In particular, for a given 
scale, it must be the same if two solutions are to describe different vacua of the same theory. 
Using the Ansatz (2), this is 

Qm X = Q+- \{k-f)F + 2Pf] . (32) 

7T 

Note that if we set L 2 B 2 = at the tip (i.e. requiring /(r = 0) = 0), then we have Q = 
Qb — Qr)3 ge m °dulo M, while the second term in (32) gives the flux contribution to the Maxwell 
charge. 



5.2 A dictionary for the charges: two puzzles and two solutions 

Our purpose is to establish using only UV data at a holographic screen whether two asymptotically- 
KS solutions describe vacua of the same theory. Any particular KS field theory is defined at a 
scale A c through a gauge group SU(Ni) x SU(N 2 ) and the associated gauge couplings (gi,g 2 )- 
The UV data of the supergravity theory consists of Qd$ x (= M), Qd3 X , j S 2 B 2l and the "stan- 
dard lore" dictionary between the supergravity UV data and the field theory is 



Ni 

4tt 2 4tt 2 



4tt; 



9i 92 
— 2~ 9s e 

92 J 



Q%r + q 

Q 



'-D3 
M ax 
D3 > 



D5 J 



-1 -* 



Bo 



7T 



mod (2 7r) 



(33) 
(34) 

(35) 
(36) 



as reviewed in [20]. We can also trade the integral of B 2 for Q^ e using 



B 2 



Max 
D3 



QPa.ge\ I n Max 
D3 )/^D5 



^ v,- 



1f/Qm x 



+ 



Bo 



^0 



(37) 



As we will see shortly, this dictionary is in fact more involved. 

All this data is defined in the supergravity solution at some UV cut-off r c related to the field 
theory scale A c . To obtain this relation, we change to a radial coordinate r such that the metric 
on the transverse six- dimensional space asymptotes to a warped conical metric: 



r/s 2 



h-^dslz + h^dsl 



(3* 
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with 



dsg ~ dr 2 + r 2 ds^i,! , r >> 1 . 
For any KS background (4), this r coordinate is related to the deformed-conifold r coordinate 



via 



^ = if 3 4' 3 ^. (39) 

The field theory cut-off A c should then be identified with the holographic cut-off r c . Note that 
from the point of view of the r coordinate, the parameter e only enters the function A from the 
Ansatz, and changing it corresponds to a rescaling of the four-dimensional metric (see (4)). 

We now run into the first puzzle, which can be expressed on the supergravity side alone. 
According to the dictionary above, since the field theory gauge group ranks depend only on 
Qd3 X but not on or one can see from equation (29) that the duals to solutions with 

different and qf but the same Q^ x have the same charges and should be dual to the same 
field theory. This is achieved by shortening the domain of integration in (31), which lowers qf, 
and by increasing qj, to compensate this. Hence, the only UV holographic data that will be 
different between, say, a solution with no mobile branes and a solution with one mobile brane 
will be the integral of B 2 on the S 2 . However, this difference is not gauge-invariant, and if one 
does not impose any infrared boundary condition on B 2 , we can see from (37) that this value 
is arbitrary, and hence nothing in the UV will distinguish between a solution with one mobile 
brane and one with no mobile brane; we expect this to be incorrect. 

One way to remedy this is to impose an IR boundary condition, namely that the integral 
of B2 on the shrunken S 2 at the tip be gauge-equivalent to zero. If so, then two solutions with 
different numbers of mobile branes and different qf will have different B fields in the UV, and 
will correspond to different theories. The only situation when the UV fields will be the same 
is when the number of mobile branes differs by multiples of M, when indeed we expect these 
solutions to correspond to different vacua of the same theory [13]. In the next subsection we will 
illustrate this in detail using our perturbation theory machinery. 

The second quandary has to do with the field theory interpretation of two solutions that have 
the same Q^ x but different numbers of mobile branes. If one is to take a holographic screen 
at r c and use the dictionary (33,34,35,36), a solution with p < M mobile branes and one with 
none will be dual to two field theories that have the same ranks of the gauge group at the same 
cutoff, but differ only in the coupling constant. Furthermore, a solution that has Q^ x — M+l 
at a holographic screen at r c will have Q^ x = M at a holographic screen placed further down 
in the infrared; this would appear to imply that a theory with rank SU(2M + 1) x SU(M + 1) 
at some energy flows at lower energies to a theory with rank SU(2M) x SU (M), then SU (2M — 
1) x SU(M — 1), which is definitely incorrect. 

A partial solution to this puzzle is given by a comment in [20], where it was noted that one 
cannot relate the UV supergravity data to field theory data at an arbitrary UV holographic 
screen. The dictionary (33,34,35,36) can only be used at special values of r c , given by the 
requirement that from the infrared up to that scale the number of duality cascades is an integer, 
or alternatively, that the value of g/ is a multiple of M. This is a stronger requirement than 
demanding that the ranks of the putative dual gauge groups are integer-valued. We will call for 
convenience the holographic screens at which one can define the dictionary "K-screens." 

However, this cannot be the whole story. As we can see from equation (37), this restriction 
alongside the requirement that B 2 be zero at the tip imply that the value of the B 2 integral at 
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the K-screen is a multiple of M, and hence the two field theory coupling constants will have 
the same values at any K-screen. Thus, at those screens (which are the only places where the 
field theory has an approximate Lagrangian description), the right-hand side of equation (36) is 
always equal to ir, and the coupling constant of one of the gauge group always becomes infinite. 
Conversely, out of the set of possible field theory data defined at a scale A c via the 4 parameters 
Ni,N 2 ,gi and g 2 , the KS supergravity solutions would only describe field theories that belong 
to a codimension-one subspace, and hence not the most generic field theory. 

In order to avoid the above-mentioned problems, equations (35) and (36) should be used 
to obtain the values of the coupling constants as a function of the corresponding energy A c . 
However, the ranks of the gauge groups given in equations (33), (34) must be read from the 
K-screen right above it. Those equations then provide the ranks of the gauge groups both at 
the scale corresponding to r c and at the scale corresponding to the K-screen above. The ranks 
do not change when one changes the position of the holographic screen by decreasing r c , unless 
one crosses another K-screen, which corresponds to a Seiberg duality in the dual theory. 

One can also ask how can a holographist tell, using purely UV data, where the K-screen lies. 
The answer is given by (36) - the screen is at the location above r c where the B 2 integral is 
gauge equivalent to zero. Hence, if the B 2 integral at the tip is zero, this dictionary gives a way 
to relate all 4 parameters of the field theory to the four parameters of the supergravity solution, 
using UV data alone. 



5.3 Baryonic and Mesonic Branches - a Perturbation-Theory Anal- 
ysis 

When the ranks of the two gauge groups are 

Nx = (k + l)M, N 2 = kM, keZ (40) 

the theory has two classically disconnected supersymmetric moduli spaces, the baryonic and 
mesonic branches [13]. For more general (Ni,N 2 ) the mesonic branch is supersymmetric while 
the baryonic branch is lifted. It is instructive to use the dictionary above together with the 
infrared boundary condition for B 2 to demonstrate in the supergravity perturbation theory 
framework we have developed that when they exist, both the baryonic and mesonic branches 
are indeed different vacua of the same theory. 

As mentioned in section (5.1), if one imposes j s2 B 2 = modulo M at the tip, then the 
function / shoud go to zero at the origin. On the other hand, we have from (150) in Appendix C 
that 

Mr = 0) = f(r = 0) + ^ = Y 7 IR , (41) 

where we have set Y§ R = since this mode diverges as 1/r 3 , and we have used the relation 
between P and M from (7). This implies that in our perturbation theory 

Qo7 = Q = - M Y 7 IR . (42) 

7T 

Setting this equal to an integer multiple of —M, leads to 5 

Q Ra f = _ £M; (43) 

=>Y r IR = ~£. (44) 

5 In our conventions the KS background has negative D3 charge. 
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Physically this corresponds to adding £ M > mobile D3-branes smeared on the tip of the KS 
solution and for each i e Z this provides the bulk dual to the E-ih mesonic branch. Let us note 
for later use that from (153), Appendix C, we get that the warp factor at the tip is 

^-^(D^-sd) 4 -" 1 - <45) 

To compare the Maxwell charges of the baryonic and mesonic branches, we must demand 
that they are defined at the same scale A c . To do so we must address the fact that the constant 
60 appearing in (39) is not gauge invariant and can be set to one by rescaling the space-time 
coordinates x^. As such one would normally fix the gauge and eliminate this constant. Indeed, 
Eq is dimensionful and just serves to fix the units which may as well be set to unity. However 
the ratio between the value of £q in two different KS vacua, such as the mesonic and baryonic 
branches, is dimensionless and physically relevant. 

This is similar to the familiar domain wall solution from one AdS vacuum to another. In either 
vacuum the AdS radius sets the units in which all other dimensionful numbers are measured but 
the ratio of the two radii is related to the ratio of central charges and is physically meaningful. 
Having said this, it is important to establish that in our Ansatz the rescaling of x^ is done by 
the constant shift in A, given in the UV by 

A = I (04 - 0i) = ~\rr + C(l/r) , (46) 

where we have preemptively used the UV boundary conditions (65) introduced below. So, 
allowing for just Y7 and Y\ to be non-zero, we can find the supergravity solution of the mesonic 
branch as a perturbation of the baryonic branch. Using (32) and (4), along with (161)-(163), we 
find that in our perturbation theory the zeroth- and first-order Maxwell charge at a particular 
radius r c » 1, is 6 

QMax = JJ^_ (r _ 1} + 8P uv + Q , T/3 s (47) 

7T 7T 

Using an expansion of e 

£ = £o(l + - + 0{Z 2 )) , (48) 
£0 

where Eq denotes that of the baryonic branch, it is apparent that if we want to stay at a fixed 
r c , then (39) requires at first order 

5r = -2 - . (49) 
£0 

Demanding that Q^ x at r c is equal for the baryonic and mesonic vacua, yields the relation 

F V UV 

Using (16) and the fact that X a = 0, we have Y 7 UV = Y 7 IR . Then, referring to (44), we have 

^ - — (51) 



'See footnote (5). 
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which is the first-order approximation to the known result Ei = €q e e7T ^ M [13, 10]. 

Now, we can find the value of the other integration constant, Y\. Using the way that e 
enters into the PT Ansatz through A, equation (4) and the UV expansions of Section (C.2) for 
A = (04 — 0i)/3 we get 



e 10 

Combining this with (50) results in an expression for Y\ in terms of Yy- 



(52) 



Yr = ^rr. (53) 

The relations obtained in this subsection can also be used to formulate the second and the 
third criteria for distinguishing between asymptotically-KS solutions. 



6 Finding the anti-D3 brane solution 

We can now summarize the necessary ingredients for identifying the candidate supergravity 
solution describing the backreaction of anti-D3 branes. Firstly, we must eliminate unphysical 
IR singularities. For many modes this is entirely unambiguous, for other modes this can be 
somewhat subtle and as such we will discuss each mode as it arises. Secondly, we demand that 
the UV asymptotics are the same as for the original KS solution which we are perturbing around. 

In total, we have sixteen integration constants but the seven physical modes (dual to seven 
gauge invariant operators) account for just fourteen of these. In addition, one is accounted for 
by the zero energy condition (9), which we use to eliminate X 5 : 

X 5 = ^ (x 2 - 2 -X^\-*-X 7 . (54) 

The zero-energy condition is necessary to completely fix the reparameterization invariance of the 
radial coordinate (see [21] for a very explicit description of this). The final mode corresponds to 
the rescaling of and for reasons discussed above this is an important physical constant which 
is given again by (52). It was pointed out in the revised version of [10] that the two vacua of 
the Klebanov-Strassler theory necessarily have different values of e. With our technology we are 
able to in fact compute the precise ratio of e in the two different vacua. 

The reader who is more interested in the end-process and in seeing or using our solution than 
in the boundary conditions we imposed to pick it out of the full parameter space of first-order 
deformations around the Klebanov-Strassler background can directly proceed to Section 6.3. 



6.1 IR boundary conditions 

We impose that the divergences in the IR for all the fields are zero, except for 4 and , the 
warp factor and 5-form flux along the brane, which should go respectively like 1/r and 1/t 2 due 
to the anti-D3-brane sources. The latter means that 05 should go to a constant. 

From the divergent term in 0§ appearing in equation (146) of Appendix C, one finds the first 
relation among Vs and Vs parameters that must be enforced: 

X 8 = l {h Q X 1 -ZPX fi ) . (55) 
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From the divergent terms in 2 we get upon using (54) that 



Y 2 IR = 0, X 6 = hoX \~ 3X \ (56) 

6 r 

Out of the divergent terms in 3 we set (after using (54) and (56)) 

Yj R = 0, X 4 = ^ X 1 . (57) 

Note that the log r/r term is automatically zero once we take into account (56). Finally, the 
divergent term in 06 requires 

Y 6 IR = . (58) 

Likewise, the other piece is zero upon using (56), (57). 

In summary, out of requiring IR regularity in all fields apart from the warp factor, we have 
obtained the following relations 

Y iR = Y iR = Y iR = Qi x 4 = ^ X 1 , X 6 = -^X 1; X 8 = h X 1 . (59) 

They are part of the relations that pick out of the full space of first order KS deformations the 
candidate solution describing the dual to a metastable state, taking into account the backreaction 
of anti-branes onto the zeroth order background. Let us move on and impose the remaining IR 
boundary conditions. 

We will now impose that there are N anti-D3 sources at the tip. The IR regularity condi- 
tions (59) yields 

5 (O) = Y 7 IR , (60) 

as in the supersymmetric case described in Section 5.3, equation (41). We require Q = N 
(cf. footnote 5), which results in 

Y 7 IR = ^p*, (61) 
where we have used (42) and (7). On the other hand, the warp factor is such that 

rto=8(?) } (w-^). (62) 

It ensues from requiring this exhibits the expected behavior for regular 3-branes (given in (45)) 
that 

X^^N. (63) 

Before moving on to discussing UV boundary conditions in the subsequent section, we note 
that inserting (63) in (10) leads to the following expression for the force exerted on a D3-brane 
probing this backreacted supersymmetry-breaking solution: 

8tt 2 2 / 3 N 
D3 ~ h 2 Q (sinh2r-2r) 2 /3 ' ( } 

This is precisely equal to the force on a probe anti-D3 brane exerted by N D3-branes that is 
computed in KKLMMT [8]. This provides further support that our IR boundary conditions are 
the right ones for anti-D3 branes. 
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6.2 UV boundary conditions 

As part of our UV boundary conditions, we impose the absence of non-normalizable modes (we 
will come back to discussing this point in section 7.2). Requiring no divergent terms in 03, 04 
as well as 05, 06 and 07 implies 

Y? v = , X 3 = , Yf v = . (65) 
Requiring no e~ T / 3 ~ 1/r terms in 02, and using (65) then determines 

X 7 = 0, X 2 = -\h X 1 . (66) 

Besides, we do not want to turn on the non-normalizable mode that shifts the dilaton, which 
would correspond in the gauge theory to changing the sum of the coupling constants for the 
gauge group. Hence, we must enforce that 

Y 8 UV = . (67) 

From (65) and (67), we see that the Maxwell charge in the UV is the same as in Section 5.3, 
equation (47). We should demand that at a given bulk radial slice r, this is the same as the 
Maxwell charge for the supersymmetric vacuum, which is in the (first) mesonic branch and has 
M — N = 4 P — N D3-branes at the bottom. Keeping in mind that e is allowed to differ in the 
two vacua, which using (39) implies that the Maxwell charges have to be evaluated at different 
r, we require that 7 

Qm x = -—(ro + 5r ms -l) + —Y 7 uv (68) 

TT TT 

8 P 2 

= (to + St! - 1) -4P + N. (69) 

TT 

Here 8t\ corresponds to the cut-off associated to the first mesonic branch. It is given by 

^i = -^, (70) 
where we have used 8 (49) and (51) for 1=1. We therefore have 



16P 2 e™ 8P x 

- 7 



TT Eq TT 

Using (52) to relate the change in e to Y^ v leads to 



ms + — Y 7 UV = N . (71) 



Y uv + _y 7 uv -N = 0. (72) 

TT 5 TT 

Note that if Y 7 UV were equal to Y 7 IR , the latter being given in (61), it would ensue that Y^ v = 
and no change in e would be necessary. However, consequent on inserting all our boundary 
conditions apart from the one associated to Y\ in (16), one finds 

8£ yUV = 8£ ylR = g ^g R (73) 

TT TT 

The shift in e can be tuned to cancel the difference in the first-order Maxwell charge Q Max 
between the anti-D3 and the supersymmetric solution. 



7 See Figure 4 below. 



8 Recall that P = \ M a! ' . For convenience we have fixed a' = 1 throughout. 
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6.3 The perturbative solution for anti-D3 branes in KS 

In summary, from the IR and the UV boundary conditions, all the integration constants turn out 
to be expressed in terms of the number N of anti-D3's at the tip of the throat. As a reminder, 
ho = h(r = 0) denotes the zeroth order warp factor of the Klebanov-Strassler solution (5) 
evaluated at the tip. Below we collect the outcome of the analysis from the previous two 
subsections: 



6 n f 1 



X 3 = , Y 3 IR = , Y? v = Z"" N . 



8.40238 
J 5 * 



X 4 =-^-iV, 1^ = 0, (74) 

i no 



7T .., , 0.70511 

6Ph 



■^N y«-0 y ^__i08244 
8P/10 ' 6 ' 6 ~ P 



*r = 0, *r'*=8^ 2^21552 

8 8/i 8 8 P 2 

All the constants in the leftmost and middle columns, with the exception of Y^, have been 
obtained by directly imposing boundary conditions in either the IR or UV. From there on, Y^ v 
was obtained from Y^ v via (72). Finally, the rightmost column was derived from the numerical 
integration which is tabulated in (16). We have not computed the value of Y[ R as it is more 
involved than the others and we do not need it, but in principle it can be done through numerical 
integration of the analytic solution (138). 

It is interesting to observe the profile of the first-order perturbation to the Maxwell D3 charge 
Qd3 X > gi ven i n Figure 3 for N = 1 (see footnote 5). Note that it does not increase monotonically. 

On Figure 4 we have plotted the total Maxwell D3 charge (i.e. the zeroth- plus first-order 
contributions) for the anti-D3-brane solution, alongside the Maxwell charge of the supersymmet- 
ric vacuum (69), the latter belonging to the first mesonic branch. For the purpose of illustrating 
equations (68)- (69), we also plot the "would-be supersymmetric vacuum" in the baryonic branch, 
that we use as a reference to measure the difference in UV cut-off, St. This branch obviously 
does not exist for N ^ 0, but it is instructive to use it as yardstick. 
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Figure 3: The profile of the first-order Maxwell charge for the anti-D3 solution, setting N = 1. 



Qd 3 (t) 

Sn <5r ms 




r 



Figure 4: Total Maxwell charge for the anti-D3 solution (blue), for the supersymmetric vac- 
uum from the first mesonic branch (red) and for the "would-be supersymmetric vacuum in the 
baryonic branch" (black dashed line), fixing iV = l,M = 3(P = |) . 



6.4 Asymptotics of the solution 

The Green's function for the KS background diverges in the IR (144), and we denote the constant 
in its series expansion around r = as jo, Eq. (145). The IR and UV series expansions of the 
solution in terms of ho, jo and X\ = N follows. 

6.4.1 Behavior in the infrared 

In the IR the solution behaves as 
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2/3 



= 33.1634 P 2 *!-^ 
= -128 (|)V^r + f (|) 5 {ho + 16P 2 j )X 1 r l + 0(7*), 

= 612.592 P 2 X X - ^ 5 P 2 Xl r + ^ * (7 ^ + 352 P 2 Jo ) X x r 2 + G(r 3 ) , 



= ^ 2 PX x -4 /iqPXx^ + O^ 3 ), 

= i/igpjJfx-y /i PX 1 + ^ ^-160/io Jo P + 10451.6P 3 ) ^ r 



3 \3 



9 \ 3 



2\ ^ 



, , //o/'.V, v 83.769 P 3 .V; r 2 + C(r d ) , 
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6.4.2 UV behavior of the solution 

As for the ultra-violet behavior of the solution, it is described by the following UV series expan- 
sions: 



08 = ~y 2 1 / 3 e" 4 ^ 3 ho X 1 (r - 1) - 288 2 2 / 3 e" 8 ^ 3 P 2 X 1 + O( e - 10T / 3 ) , 

2 = -418.571 e~ T P 2 X l + — 2 1/3 e~ 7r/3 h X 1 (1 + 8r) + C(e" 3r ) , 

3 = - f 2 1 / 3 e" 4 ^ 3 fco X x + 2 e" 2r (1186.08 - 418.571 r) P 2 X x - ^ 2 2 / 3 e" 8 ^ 3 P 2 X x 
+ O(e- 10r / 3 ) , 

0i = 428.85 P 2 Xx + ^ 2 1/3 e~ 4T/3 h X x - e~ 2r (1325.73 - 837.143 r) P 2 X x 
9zL 

+ ^ 2 2 / 3 e" 8 ^ 3 P 2 (29 + 40 r) X 1 + 0(e^ 3 ) , 
5 

5 = 312.743 P 3 Xi + e- T (-1361.84 + 418.571 r) P 3 X x - 4 2 1/3 e" 4r/3 h P X x (l + 8r) 

+ 2e" 2r (1361.84- 837.143 r) P 3 X 1 + C(e" 7r/3 ) , 

6 = 312.743 P 3 Xx + e~ T (1361.84 - 418.571 r) P 3 X x - 4 2 1/3 e~ 4r/3 h PIi(l + 8 r) 

+ 2 e~ 2T (1361.84 - 837.143 r) P 3 X a + 0(e~ 7r/3 ) , 

7 = e" r (943.269 - 418.571 r) P 3 X x 

- ^ 2 1/3 e" 7T/3 /io P (1199 + 80 r (1 + 10 r)) X x + £>( e - llr/3 ) , 

<A - 171 54P 3 X + 4 2 1 / 3 e- 4 -/ 3 /. (7 + 32r)X 1 625.486P 2 X 1 

4 - 171.54P X X + 3(4^31) (4r-l) + ° (e } ' 

7 Additional Comments 



83) 
84) 

85) 

86) 
87) 



89) 
90) 



Having solved for the full space of linearized perturbations around the Klebanov-Strassler back- 
ground, we now discuss other solutions that we easily obtain as a by-product of our analysis, as 
well as other possible interpretations of our results. 

7.1 Relation to previous works 

The first attempt to construct the a linearized antibrane solution in the UV region alone was 
[14], which studied several of the SU(2) x ST/ (2) x Z 2 -invariant modes around the Klebanov- 
Tseytlin (KT) background [22]. Since the KT solution is a subset of the parametrization (l)-(2) 
given by 

y(r) = 0, fc(r) = /(r), F(r) = P , (91) 

in our setup we can understand the perturbations around KTas solutions of a reduced system of 
first-order differential equations in the Borokhov-Gubser formalism. The details of this analysis, 
as well as the relation with the notations of [14] can be found in Appendix D. The "backreacted" 
KT solution contains some integration constants that cannot be fixed by infrared boundary 
conditions, and hence we cannot relate them to the constant X%, which is proportional to N . 
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We can directly compare the UV expansion of our full KS solution (83)- (90) to the perturbed 
KT solution of [14] and we find the following crucial discrepancy: The correct UV expansion 
has terms of order 0(r~ 3 ) in (84,87,88,89) while the first non-trivial terms in the solution of [14] 
are at C(r~ 4 ). 

In hindsight this is not so surprising, since [14] only considered a subset of the modes, and 
furthermore, the KT solution precisely agrees with the UV limit of the KS solution only at 
leading order. At subleading order the KT solution has an ambiguity which can be fixed to 
agree with the UV limit of the KS solution but then the lower-order perturbation theory around 
each solution quantitatively differs. For this reason, we conclude that one cannot derive the 
correct UV expansion for the anti-brane solution by starting with the KT geometry. Another 
problematic issue with the Ansatz made in [14] is that, as we have explicitly demonstrated in 
this work, the anti-D3-branes turns on modes which are outside of the truncation, so it is not 
consistent to restrict oneself to this subset of mode. We refer the reader to Appendix D for a 
thorough analysis of those issues. 

7.2 Gaugino masses 

As an additional outcome of our analysis, we can easily identify other interesting solutions that 
correspond to different deformations of the dual gauge theory. In particular, we can construct 
a solution in which the non-normalizable UV modes X 2 and X7 are turned on. They decay as 
1/r, and are associated to operators of dimension A = 3, which correspond to deformations 
by gaugino mass terms for each of the gauge groups, Tr(AiAi ± A 2 A 2 ). We will identify a one- 
parameter subfamily for which Q^ x approaches the same constant value in the IR and in the 
UV, and therefore for which the parameter e does not need to be modified. 

The boundary conditions we have to impose are exactly the same as before, except that now 
we do not require (66). Relaxing these, we find that the leading terms in the IR expansions 
are not modified, and the value of 05 at the origin is still given by (60), together with the 
relations (61)), (63) 

5 (O) = Y 7 IR = (92) 

By using the UV/IR relation (16) we get that in the UV 

5 (oo) = Y 7 UV = 154.299 P 3 X x - 19.5477 P(2X 2 + PX 7 ) . (93) 

Imposing 5 (0) = 05 (00), we thus see that for the family of solutions 

2X 2 + PX 7 = 5.05767 P 2 X x (94) 

we get that the first order Maxwell D3 charge at infinity is the same as that of the supersymmetric 
vacuum with the same e as for the original KS background. The profile of the perturbation to 
the D3-brane Maxwell charge is shown in Figure 5, where it is plotted as a function of N using 
the condition from equation (94). 

We also note that by setting X\ = 0, i.e. requiring that no anti-D3 brane be present at 
the origin, we obtain a family of solutions parametrized by the constants X 2 and X7 which in 
the dual gauge theory describe soft supersymmetry breaking due to gaugino mass terms. This 
solution encompasses the one built in [23], which corresponds 9 to the family X 2 = PX7 . 

9 The constant X in [23] is then related to X 7 by X — — \X 7 and their parameter fi is such that [i = 48 2 1 / 3 PX 7 . 
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Figure 5: The profile of the first order Maxwell charge for the solution with gaugino masses turned 
on, satisfying the constraint (94) (blue solid line). The plot is for N — 1 and X 7 = 1/(24 2 1 / 3 P 3 ). 
The red dashed curve is the profile for N — 0. 



7.3 Other UV boundary conditions 

In section 6 we have identified the anti-D3 backreacted solution using one of the three criteria 
to distinguish asymptotically-KS supersymmetric solutions that we have put forth in section 5. 
The resulting solution has a different scale parameter Yi than its supersymmetric counterpart, 
and if the criterion that the NSNS B2 field be zero at the KS tip is the correct one, then, 
putting aside concerns about the subleading singularity and about backreaction, the anti-D3 
perturbative solution we have constructed describes a metastable stable of a supersymmetric KS 
field theory, and would be the first metastable solution constructed in supergravity. 

However, we can also ask whether this result holds if one imposes the other criteria, or if one 
insists, perhaps with a view towards embedding the KS solution in a compact setting, that the 
UV scale parameter Y\ be the same as in the supersymmetric theory. It is not hard to see that if 
one imposes the criterion that the H% integral only jumps by integer units, one finds again that 
Y"]_ has to change; the anti-D3 solution is identical to the one we have written down above, 
and would be dual also to a metastable field theory vacuum. 

If one on the other hand imposes the criterion that two vacua of the same theory must have 
a Y^ v related to Yj JV as in equation (53) (which also distinguishes between various supersym- 
metric KS vacua), or imposes the requirement that the UV scale must be the same as in the 
supersymmetric theory, then the resulting solution will have a different IR Maxwell charge than 
the one inferred from the UV data (essentially because antibranes give rise to negative charge 
dissolved in flux in their vicinity, as shown in Figure 4, and if one cannot make the throat longer 
to compensate for this, this charge will be visible at infinity). As a result, the relation between 
the force on a probe D3 brane and the anti-D3 charge of the background will not be the one of 
[8] . If one then insists that this relation does not receive corrections at first order in the number 
of antibranes, as suggested by the no-screening results of [7], then the anti-D3 solution must 
have a nontrivial 1/r mode turned on, of the type presented in the previous subsection, such 
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that the contribution to the charge dissolved in flux from the antibranes is canceled by the con- 
tribution from the X 2 and X 7 modes. The value of this non-normalizable relevant perturbation 
can be easily read off from our analysis. Interestingly enough, such modes were argued in [24] 
to be present when a KS solution is embedded in a stabilized flux compactification, and it is 
interesting to see if the relation between the anti-D3 charge and the strength of this mode that 
one can find here has any relevance to this analysis. 
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A First order deformation around a supersymmetric back- 
ground 

A.l Our approach: from second order to twice as many first order 
differential equations 

The method introduced by Borokhov and Gubser [4] to find the set of first order perturbations on 
top of a supersymmetric solution depending on a single radial variable r, relies on the existence 
of a superpotential W whose square gives the potential obtained by reducing a supergravity 
Ansatz: 

1 , dW dW 

v ^S G "arW (95) 

The fields a (a = l,...,n) are expanded around their respective supersymmetric background 
values 4>q, 

r = <Po + <Pl(Z) + 0(Z 2 ), (96) 

where Za = {X a , Y a ) represents the set of perturbation parameters and <P1 is linear in them. 
The method amounts to splitting n second-order equations into 2 n first-order ones, out of which 
n of them (those for the conjugate momenta £ a ) form a closed set. The defining equations for 
the modes £ a are 

Ha = GMo) ( #i - M b d (<f> ) ti) , M\ = \4rA G bc SI • (97) 



dr *™-rij , a 2d<P d \ Oft 

They measure the deviation from the BPS flow equations, i.e. they are non-vanishing only for 
non-supersymmetric solutions. The set (£ , <p a ) satisfies the equations: 

^+£ 6 M*> a (0 o ) = 0, (98) 
' h) " ] M\{<P,)<P\ = G ab £ b , (99) 



dr 
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where (99) is simply a rewriting of (97) whereas the equations in (98) imply the equations of 
motion [4]. Additionally, the functions 6 should obey the zero-energy condition 

6 ^ = o , (100) 

which stems from gauge-fixing the additional degree of freedom corresponding to reparametri- 
sations of the radial variable (as explained very clearly in [21]). 

The n integration constants arising upon solving (98) are branded X a , while those associated 
to (99) are identified as Y a . 



A. 2 First-order equations for the deformations around KS 



Let us now review how the Borokhov-Gubser method is implemented for studying perturbations 
around the Klebanov-Strassler solution. There are eight functions in the Papadopoulos-Tseytlin 
Ansatz written in (l)-(2), a = (x,y,p,A,f,k,F,Q), which is a consistent supersymmetric 
truncation of type IIB [25, 26]. Their zeroth-order values are available above in (4). 

The field-space metric entering equation (97) is computed out of the kinetic terms arising 
from the IIB reduction 



G, 



ab ' 



L'a lib 



12 i 1 12 , a '2 a A>2 , ^ ^/2 

x + -y +6p —6A + - <p 
2 4 



1 

+ 4 C 



-$>-2x 



( e -2y f ,2 + e 2y k ,2 + 2e 2Z F ,2j 



(101) 



The superpotential is found from the corresponding potential appearing from the reduction of 
the PT Ansatz: 

W(0) = e 4A ~ 2p - 2x + e 4A+4p coshy +^e 4A+4p ~ 2x (/ (2 P - F) + k F) . (102) 

In order to solve the system of equations (98), (99) for the modes 6 and 0", we find it convenient 
to rotate to a different basis (6, a ), defined as follows in terms of the original fields: 

6 = (36 - 6 + 6, 6, -3^ + 26 - 6, -36 + 6 - 26, 6 + 6, 6 - 6, 6, 6) 

4> a = (x-2p-5A, y, x + 3p, x-2p-2A, f, k, F, $) . 
In the order we solve them, the system of first-order equations for the 6 (98) reads 
g = e -a*o [2P/ -Fo(/o-*o)]fi, 
I4 = -e~ 2x ° [2Pf -F (f -k )]i u 
g = - l -Pe- 2 *°~ 
1 



is 



6, 

-2x n 



(103) 
(104) 

(105) 
(106) 

(107) 



-6 - o e (P-F ) 6 



1 



- sinh(2 y ) 6 - cosh(2 y ) 6 + - e~ 2x ° (f - k ) 6 , 

6 

(Pe 2yo - sinh(2y ) F ) 6 + {Pe 2yo - cosh(2y ) F ) 6 + \ (fo ~ h) 6 



3e 



-2x -6p 



6+ [5. 



-2x -6po 



-2x 



(2P/o-F (/o -**))] 6 



6 cosh ?/o + g sinh|/o (26+6 + 6 



+2 



( Pe 2yo _ cosh(2y Q ) F ) 6 + (Pe 2s » - sinh(2y ) F ) 6 



(108) 
(109) 

(HO) 
(HI) 

(112) 
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Particularized to the deformation around KS, the system of equations is 



_ 4e -4^ -4po<£ 8; 

-coshyo02-2e- 4Ao - 4po 6 
—3 e ~ 6po_2:Eo 



3 - sinh y 02 - e 
o 



-4A0-4P0 



96 + 56 + 2^4 , 



2e 



-6po-2a;o 



i 3 -sinhyo02 + ge- 4A °- 4po (Ii + 3| 4 j • 



e 22/o ( Fo _ 2 p) (2 02 + 8 ) + e 2yo fa -2e- AA °- 4po+2xo+2yo (| 5 + f, 



-2»o 



F n 2 



2 e -4A -4 P0+2 0=0-2^0 ^ 5 -^ 6/ . 
' -2x \ f /n r> T?\\UT?~\(nl n 2 el \ i ^ „-2a; 



ge^I/o (2P-F ) + A; Fo] (^2 X - 2 3 - 5 4 J + ^ e^ (2P- 



1 

+ 2 e ~ 



2 "F a ^ + -e- 2 " (k a -h) & 



-4A -4p ^ 



(113) 
(114) 

(115) 
(116) 
(117) 
(118) 
(119) 

(120) 



B The analytic solution space of deformations around KS 

Here we provide for handiness the solutions for the £ a 's and </> a 's, in the order in which they were 
solved in our previous work [6]. Of main interest are the <p a but we first had to solve for their 
"conjugate momenta" £ a sourcing their equations. 
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B.l Analytic expressions for the £ a modes 



6 

is 
u 

ir 

is 
6 



X, h(r) , 
5 X x h{r) 



32 



3 
128 



P 2 X 1 csch 2 r (sinh r cosh r — r 



, 4 / 3 



P 2 Xi (sinhr coshr — r) j(r) + 2X3 (coshr sinhr — r) 



9 

-Xx /i(r) + X 4 , 
- 1 ^X 1 j(r)+X 5 , 



Ae(r) 



cosh r sinh r — r 



sinhr 

coshr , . , 
— r^A 6 (r) + 
smh r 



2 sinhr A7(T) ' 
-3 + cosh 2 r + 2 r coth r 



P (r coth r - 1) cothr£ 5 - P 



4 sinhr 

r cothr 



A 7 



sinh r 



^ 8 



6 



X!/i(r)+X 8 



2 1 , 

— X 3 r cosh r H — X 4 cosh r + P X 6 cschr (coth r — r 

3 3 

+P X 5 cschr (l — 2 r cothr + r 2 csch 2 r) + X 2 sinh r 
+ ^PX 7 (-2r coth 3 r + csch 2 r + r 2 csch 4 r) sinhr 



csch 2 r) 
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X 1 



3 csch 3 r h(r) (6 r — 5 sinh 2 r + sinh 4 r) 



(121) 

(122) 
(123) 
(124) 

(125) 

(126) 

(127) 



+2P 2 csch 5 r ( - 15 + 24 r 2 + 16 cosh 2 r - cosh 4 r - 32 r sinh 2 r + 4r sinh 4 r) 



4 sinh 2 rj(r) — 6 (coshr sinhr — r 



,V3 



(128) 



where 



1 IX 

A 6 (r) = X 6 + - (-r + cothr - r coth 2 r) £ 6 (r) + - /i(r) , 
z 6 P 

16 

A 7 (r) = X 7 — csch 2 r £ 5 (r) + — P X 1 csch 2 r (coshr sinhr — r 

o 

64 

+ -PX l3 (r). 



1/3 



(129) 



(130) 



B.2 Analytic solutions for the 0f's 



Holding our breath, we recap the analytic solutions for all eight </>f modes found in [6] : 
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r J sinn u (cosh-u sinhu- u) 1 

2 ~ 16 64 

+ p Ht) i 5 (r) + — X 1 cschV (cosh r sinhr - r) 1/3 /i(r) + — X 1 h(r) j(r) 

64 /" r ( sinh 2 tt + 1 — u coth it) „ s 

+ y Xl / - U r h h — w* h{u)du > (131) 

<J J sinn it (cosn-u smhu- u) ' 
2 = cschrA 2 (r), (132) 

* = 53^7=27 ^ < 133 > 

4>i=Y± + — X 4 j(r) - - 03 (r) — X 3 / (coshw sinh w - w) 1/3 dw 

+ coth -u A' 2 (-u) du — coth rA 2 (r) H — ^^P 2 Xx ^ csch 2 w (cosh it sinh it — -u) 2 ^ 3 du 



H — — P Xi / (cosh u sinh it — it) ' j(u)du— — Xi / jj^du, 

81 7 27 J (cosh it sinh it — it) ' 



(134) 



$ 5 = i sech 2 (r/2) [r + 2 r cosh r - (2 + cosh r ) sinh r] A 5 (r ) + - - ^—^ ^ A 6 (r ) + A 7 (r) , 



(135) 

r / 1 \ 1 

06 



r ( 2 - l-coshr 1 " t " i ''"' r/2) 



A ^( r ) + 1 iT" A e( T ) + Mr) , (136) 

1 — cosh r 



07 = (— coshr + t cschr) A 5 (r) — cschr A 6 (r) , (137) 

t 1 Ur 16^ r M«) 2 , „ [ T ( u cothu- 1) csch 2 wA 6 (w) , 

04 = -rr^ \ Y A - — X l / ^ — du + 32P / ^ '- -f^ dw 

'H r ) 3 J (cosh it sinh u — u) 1 J (cosh u sinh it — u) 1 

^du+—P / (it coth u — 1) csch 2 w (cosh u sinh it — -u) 1 ^ 3 

(cosh it sinh it — it) 5 y 

5A 5 (iz) + 2P (-0 1 (iz) + 3 H)] du}, (138) 



x 



where 
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/ T (—2 u coth 3 u + csch 2 M + u 2 csch 4 -u) sinh 2 u 
~ wT^ du 
(coshM sinh-u — u) 1 



nn „ „ r coth u — m csch 2 M , „ „ [ T 1 — 2 m coth u + u 2 csch 2 M , 

-32PX 6 / — du-32PX 5 / 27^ M 

J (cosh u sinh m — u) J (cosh m sinh u — m) 7 

32 /" r cosh-u sinh-u 64 /" T u cos!im sinhw 

- — X 4 / du + —X 3 / ^ gJm 

y (cosh -u sinh -u — -u) ' <J J (cosh -u sinh m — -u) 7 

„ . , . , ,i/Q 8„^ [ T 6 m — 5 sinh 2m + sinh 4 -u , . . , 
-48X 2 (coshr smhr-r) 1/3 + -X 1 / - — - ^h(u)du 



9 J sinh 2 u (cosh-u sinh-u — u) 2 ^ 3 
32 2 _ f T -15 + 24m 2 + 16 cosh2w - cosh4w - 32m sinh2w + 4m sinh4w 



-—P'X 1 / : : m : du 

9 J sinh u (cosh -u sinh u — u) 

64 9 /" r -15 + 24m 2 + 16 cosh2w - cosh4w - 32m sinh2w + 4m sinh4w . 

+ 27 P Xl / ■ , 2 , r rr ^ j (m) g?m , 

z < J sinh m (coshM smhw- u) ' 

(139) 



A 3 = Y 3 



80 

y 



32 r T 112 /" r 
— ^~-^4 / (cosh m sinh m — m) 1/<3 c?m %~ ^ / ( c °sh.u sinh u — u) 1 ^ 3 h(u) du 

/T pT 
(cosh m sinh u — u) 1 ^ 3 ^(m) g?m + 2r coth rA 2 (r) — 2 / m cothM A' 2 (m) g?m , (140) 

1 ^ / , ,s , 9 7 / \ ™ ^ f T (u cothM — 1) csch 2 M ~ . . , 
A 5 = Y 5 - - P (r coth r - 1) csch 2 r 8 (r) - 32 P / } ^ — g?m 

^ 7 (cosh m sinh m — m) ' 

1 f T ( T 2 + cosh 2 m 

+ - X 7 / csch 4 w [2 m (2 + cosh 2 u) — 3 sinh 2 u\ h(u) du — X 6 / -. h(u) du 

4 y J sinh m 

„ I cosh. 7~ Slllll T~ T~ 

csch 2 M [—3 cothw + M (2 + 3csch 2 M)l h(u)£ 5 (u)du P A 2 (t) 

2 sinh r 

1 f T X f T 

+ - P / csch 4 M (coshM sinliM — m) A 2 (m) g?m — ^-^ / (2 + cosh 2 m) csch 4 M h 2 (u) du 

16 f T 

+ — P Xi I csch 4 M [2 m (2 + cosh 2 m) — 3 sinh 2 m] j(u) h(u) du 
4 f T 

+ -PX 1 / csch 6 M (coshM sinh m-m) 1/3 [2 m (2 + cosh 2 m) -3 sinh 2m] h(u) du , (141) 

A 6 = y 6 _ I p [-t + cothr + t (-2 + r cothr) csch 2 r] 8 (r) 

/ T \—u + cothM + m (—2 + m cothM) csch 2 Ml ~ , 
r h • h du 
(coshM smhM — u) 

1 [ T 

+ -X 7 / [cosh 2 m + csch 2 M (3 + 2 m 2 — 6 u cothM + 3 m 2 csch 2 M)] h{u)du 

+ X e J csch 2 M [3 cothM — m (2 + 3 csch 2 M)] h(u) du 

+ J [l + (3 + 2 m 2 — 6 m coth m) csch 2 M + 3 m 2 csc1i 4 m] h(u) I5 (u) du 
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P [2 coth 2 r (-1 + r cothr) + csch 2 r - r 2 csch 4 r] A 2 (r) 



P / [2 coth 2 m (—1 + u cothw) + csch 2 M — u 2 csch 4 u] A! 2 (u) du 



T f csch 4 M [—2 u (2 + cosh 2 it) + 3 sinh 2 m] 



12P 



fc(it) H Pcsch 6 M 

36 



8j(u) sinh 2 M + 6 (cosIim sinliM — m) 1//3 — 28 + 32 m 2 + (31 + 16 m 2 ) cosh 2 m 



— 4 cosh 4 u + cosh 6 u — 48 it sinh 2 m 



/i(m) g?m 



A 7 = Y 7 + P [-r + coth r + r (-2 + r coth r) cschV] </> 8 (r 
" T [—it + cotliM + u (—2 + m cothw) csc1i 2 m] ~ 



+ 64P 



&{u) du 



(coshw sinhw — m) 2 ^ 3 
+ Xj J [— 1 + (— 3 — 2 u 2 + 6 u cothw) csch 2 u — 3 u 2 csch 4 u] h(u) du 

+ Xq J csch 4 M [2u (2 + cosh 2 u) — 3 sinh 2 m] h(u) du 

[— 2 — 2csch 2 M (3 + 2-u 2 — 6m coth-u + 3m 2 csc1i 2 m)] h(u)^(u)du 
-PcschV (l - 2r cothr + r 2 cschV) A 2 (r) 
+ P J csch 2 M (l — 2 it coth m + m 2 csc!i 2 m) A' 2 (m) g?m 



csch 4 M [2 m (2 + cosh 2 m) — 3 sinh 2 u] 
6P 



h(u) 



Pcsch M 



8j(u) sinh 2 M + 6 (coshu sinhu — u) 1 ^ 3 

[-9 + 16m 2 + 8 (l +m 2 ) cosh2M + cosh4M - 24m sinh2M] J h{u) du . 



(142) 



(143) 



C IR and UV expansions of our analytic solutions 
C.l IR expansions 

The IR behavior of the modes is obtained by Taylor expanding h, j and the integrands in (131- 
138), performing the indefinite integral over r (instead of the integral from 1 to r), and adding 
an integration constant Y^ R (since the conjugate momenta £ a do not involve integrals other than 
h and j, we do not have to introduce a second set of integration constants X IR different from 
the one used in (121)-(127)). 
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The IR expansions of h and j are given by 



16 (2 



3 



hm = h - f I | 1 P 2 r 2 + 0(t 3 ) , 
where 

h = 18.2373P 2 , jo = 0.836941 . (145) 

In the order that those equations were solved and to the order of expansions that we need, 
the IR asymptotics of the (fi a modes are given by 

1 32 /2\ 5 

08 = -yfgJ (-h X 1 + 3PX 6 + 9X 8 ) + Y 8 IR + O(r), (146) 



02 = \Yl R + y^(y 3 (h X 1 - 3(X 4 + 2PX 6 )) j + 8 3 (-6X 2 + 4X 3 + 6PX 5 + 9PX 7 ) 
+ 0{t) , (147) 



~ 3Y 3 IR 1 Y* R 3Y 3 IR 4/2\3 r /2V nv 

+ — I I ffl 3 Cx>*i - 3(X 4 + 2PX 6 )) I - 0( r) . CI 18) 



r \ 3 V 3 



1 Y IR 1 / 4 /2\ s 

' * + _(_ 2 y/* + ^L---^-J (4/io*i - 3(5*4 + 12PX 6 )) 



r 3 2 



+ ^ ( " f (I) 5 {h ° Xl ~ 3(X4 + 2PXe)) ) + 

+ logr ( y 3 (-6X 2 + 4X 3 + 6PX 5 + 9PX 7 )^) + O(r) , (149) 



VlR 

k = \ + Y 7 IR 



, , PK, /i? y/ fl 1 o /2\ 3 1/ i2 9 \ 

t I 1 ^ + 36P^ Xl_4 V3j PX 4+g (-322333P 2 + / i0 )X 6 - 8 2333PX 8 

+ r 2 log r ( -- I V P (/i *i - 3(X 4 + 2PX 6 )) ] + (9(r 3 ) , (150) 
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2Y* R + 



8 



3 V 6P 



+ 



f + K/"- 2 ^- 1 ? 8 (t\ '' KPX, + + 16 (f| 1 P.Y, 



9 V 3 J 27P V 3 



log r (J J 3 P (ZioXi - 3(X 4 + 2PX 6 )) ) + 0(r) 



(151) 



PF« , F/* , 64/2\5 1 , 2 „ 8/2\ 1/3 _, 1 




+ r — + — + — - h PX 1 + /i^Xi - - - PX 4 + -/i X 6 - 16 

lo 6 9 1 3 / 54P 3 1 3 / 9 




16 /2V 
9" 1 3 



P (ZioXi - 3(X 4 + 2PX 6 )) + 0(t 2 ) , (152) 



-6PF 3 /i? - 18Y 6 IR - 27Y 7 IR + -^/igXi - 12h X 6 ) ) YI"+0{t) 



(153) 



Note that the constant term in 2 and the logarithmic term in <pi are identically vanishing once 
we impose the zero-energy condition (9). The relation between the constants (X, Y IR ) used here 
and those that first appeared in [5], which we denote (X IR ,Y IR ), is summarized in the next 
subsection. 

C.l.l Relation to the IR series expansion of [5] 

The relation between the X a , Y^ R integration constants in this paper and the IR integration 
constants in [5], which we call X^ R ,Yj R depends ho and jo, whose numeric values are given 
by (145). We have given by 

X( R = h X 1 , X IR = i- (-9/i + 16j P 2 ) X l + l -X 2 + ^X 4 (154) 

54 I o 

vIR — 32j P 2 y i 1 V Y IR — h Y 4- Y 

X IR = -m'X 1 + X 5 , X IR = |(-^ + 16 (2V33V3 + Jo ) Xl -X B - 2X 6 , 

X IR = -^fX, - \X 7 , X IR = + I (2V3 3 V3 + 2jo ) P 2 Xi _ p X5 _ ?x 7 + X 8 
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and 

Yl R = Yl R fora^6,7, 



16/ln ^ 3 



C.2 UV expansions 

The UV asymptotics of h(r) and j(r) are 

h uv = 12 2 1 / 3 P 2 (4r - l)^ 3 - — 2 1 / 3 P 2 (12 - 85r + 25r 2 ) e - 10 ^ 3 + 0( e - 16T / 3 ) 

= -o^e- 4 ^ 3 - 1 2^(3 + 10r)e- 10 ^ + G( e -^/ 3 ) . (156) 

The UV expansions for the fields a are obtained by performing an indefinite integration of the 
UV series of the integrands as in the IR case. We call Y^ v the Oth-order term in the expansion 
for the field <p a (or A a if the former is written as a product of the homogeneous solution times 

K) 

8 = Y s uv + 12 • 2 1 / 3 e- 4r/3 (p(-l + 4r)(2X 5 + X 7 ) + 8X 8 ) + 0(e- 8r / 3 ) , (157) 

2 = -8 • 2 1/3 e' r/3 (QX 2 + (6 - 4r)X 3 + 2X 4 + 9PX 7 - 6PrX 7 ) + 2 e~ T Y^ v 

+ 0(e^ /3 ) , (158) 

03 = -5 • 2 1/3 X 3 e 2r/3 - \ • 2 1/3 e- 4r/3 (l08X 2 + (336 - 137r)X 3 + 48X 4 

- 108P(-3 + r)X 7 ) + 0(e- 2r ) , (159) 

fa = Y? v - 10 • 2 1/3 X 3 e 2r/3 + - • 2 v V 4r/3 (324X 2 + (528 - 316r)X 3 + 1UX 4 

3 V 

+ 81P(7 - 4r)X 7 ) + 0(e- 2T ) , (160) 
vuv 

k = ~\^- Y 5 UV + Y 7 UV + r(2Y 5 uv - PY 8 UV ) 

+ 6 • 2 1/3 e- r/3 p(QX 2 + (21 - 4r)X 3 + 2X 4 + 21PX 7 ) 

+ l -e' T ((5 - Ar)Y 5 uv + 4Y 6 UV - 2P(-1 + 2t)(Y 2 uv - Y 8 uv )^j 

+ 12 • 2 1/3 e" 4r/3 p( - 12(1 + r)X 2 - 15X 3 - 4X 4 + 2r(X 3 + 4rX 3 - 2X 4 + 6PX 5 ) 

+ 3P(-3 + r + 4r 2 )X 7 + 6(PX 5 + X s fj + C(e~ 2r ) , (161) 
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Y Jl_ e r _ yvv + Y uv + T[2Y uv _ pyuv. 



- 6 • 2 1/3 e~ r/3 p(6X 2 + (21 - 4r)X 3 + 2X 4 + 21PX 7 
+ \e- T ((-5 + 4r)lf y - AY 6 UV + 2P(-1 + 2r)(F 2 c/y - Y 8 uv )^j 
+ 12 • 2 1/3 e" 4r/3 p( - 12(1 + r)X 2 - 15X 3 - 4X 4 + 2r(X 3 + 4rX 3 - 2X 4 + 6PX 5 ) 
+ 3P(-3 + r + 4r 2 )X 7 + 6(PX 5 + X 8 )) + C(e" 2r ) , (162) 



uv , 

e r + 18 ■ 2 1/3 e- r/3 Pf - 6X 2 + (-9 + 4r)X 3 - 2(X 4 + P(5 - 2r)X 7i 



2 



+ e~ T ({- 1 - + 2r)Y 5 uv - 2Y 6 UV + P(Y 2 UV + 2tY 2 uv - Y 8 uv )^j + 0(e~ 7T / 3 ) , (163) 

~ _ Y™ r 4W3 _ 8-2V 3 (2r + l)X 3 /3 2^ _ ^ 

^ 4 12 • 2V3(4r - 1) 4r-l + 5 P 2 

, 4-2 2 / 3 (12 - 85r + 25r 2 )^ y /3 , 2 1 / 3 /3 / 

" P(4r - 1) + 1125(4r-l)2 6 + (47=1) 6 v 18(? + 8r)X2 

+ 32(2r + 1)X 4 - 18P(7 + 8r)X 5 - 9P(23 + 8r + 32r 2 )X 7 - 72X 8 

40803 - 170884r + 161 120r 2 - 332800r 3 )X 3 \ ml 2rX ..,„,, 

H ; : — + O e~ 2r ) . (164) 

375(4r-l) / v ; K J 

C.2.1 Relation to the UV series expansion of [5] 

The relation between the X a , Y^ integration constants used in the present paper and the UV 
integration constants introduced in [5], which we denote here X^ v ,Y^ V , goes as follows: 

1 „ 



xf y 


= -2UP 2 X 4 , 


X™ 


= - 2 x 2 + 




= 2^3 , 


n v 


= x 4 , 


xr 


= Xb, 


vuv 


= -x 5 - 


X 7 UV 




vuv 


= -PX 5 



—y\-A . 
6 ' 



(165) 



and 



A 6 j 

-X 7 + X 8 (166) 



Y a UV = Y a uv fora^6, (167) 

(168) 



yuv _ y uv — Y uv + -Y uv 

6 6 2 ~ ^ 8 



D The Klebanov— Tseytlin perturbation 

In our parametrization of the metric and the fluxes (l)-(2), the Klebanov-Tseytlin background 
corresponds to the subset defined via (91). At zeroth-order the fields cf> a KT obey the flow equations 

« T _v«L i (169) 



dr 2 d<p% 



KT 
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with the superpotential 



W(<f>) = e iA ~ 2p ~ 2x + e 4A+4p (1 + Pe~ 2x f) . (170) 
These equations are solved by 

A = --log(h K T(r)) , 

l, f h KT ( r y \ 



1, / 48 ■ 2 1 / 3 

PO = 77 log 



6 \kKT( r ) r 

/ = P(l-31ogr), 
$o = 0, 

where the warp factor ft^T of the Klebanov-Tseytlin solution takes the following expression: 

12 ■2 1 ^P 2 {Yl\ogr - 1) 



(172) 



In order to match the UV asymptotic of the Klebanov-Strassler modes (4) we should use 
the relation r = e t//3 , while the perturbation in [14] corresponds to changing the origin of the 
log as follow: logr — > logr — |. This is equivalent to changing e. The relation to the functions 
a(r),b(r), k(r) used in [14] is the following 

a(r) = — ^(x(r) + log 6) , (173) 

b(r) = —3p(r) — x(r) — log 6 + log(3"\/6) , (174) 
k(r) = 6/(r) , (175) 

whereas the constant M is related to our P as 



M=-18P. (176) 

If we consider linearized deformations around (171), it is quite simple to solve analytically the 
linearized equations (98)-(99) for the five £ a and five 0° modes. In this way we get a solution 
which contains terms up to the order r~ 8 . Henceforth, as a bonus, we obtain the perturbation 
around the KT background that includes the mode responsible for the force on a probe D3 brane 
discussed (but not worked out quantitavely) in [14]. The results for the UV modes expansions 
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are as follows: is the following 



-288 • 2 2 / 3 P 2 X 1 72-2 1 / 3 PX / (l + 41ogr) 96 • 2 1 / 3 X$ 
~8 + ~4 + -4 + r * ' 

1152- 2 2 / 3 P 2 Xi , 16- 2^X4 F 3 , . 
^ ~ ~ 10 ■ 2l/V *3 ^ + ^ , (1T8) 



24-2 2 / 3 P 2 X 1 (29 + 1201ogr) _ 2Q + 4 ■ 2^ 3 X 4 _ 2^ + ^ 

7 2. 2 2/3p3 Xi ( 121ogr _ !) 144-2 1 / 3 P 2 X / (l + 31ogr) 



(179) 



79 ol/3 p y 

+ r4 * - 3PF 8 log r + Y) , (180) 

48-2 2 / 3 P 2 Xi(-67- 72 logr +2880 log 2 r) 18 ■ 2 1 / 3 PX f (ll + 241ogr) 

94 — 



25r 8 (121ogr - 1) r 4 (121ogr-l) 
72-2V 3 X* 16 ■ 2 1 / 3 r 2 X 3 (l + 61ogr) 2 • 2V 3 X 4 (24 logr - 5) 
+ r 4 (l - 12 logr) 12 log r - 1 + r 4 (121ogr - 1) 

2^ 217 + Pr 4 F 4 8r 3 (301ogr-7) F„(3 + 121ogr) 
5 P(121ogr-l) 75r 6 (121ogr - 1) 2-241ogr ' [ ' 

where 01,03,04 are defined as in (8), and 0/ and 0$ are respectively the perturbations to the 
function / (=k) and the dilaton. 

In order to compare to the full KS solution (83)-(90) we should identify 

Xf = X 5 , X$ = X 8 — PX 5 , (182) 

and indeed after replacing X 4 , X 5 and X 8 with the boundary conditions given in (74) the r -4 
and r~ 8 terms agree. 

By rescaling the radial coordinate r we can compare to [14]. We get the following relation 
between their parameters S, 0, and our Xf, X$ 

S = -96 2 l/3 PX f , = 24 2 1/3 (7PX / + 4X$) . (183) 

Note that the IR boundary conditions relate Xf and X$ to N 

X f = l — ttN , X^ = ^—TiN. (184) 

As a result, those parameters cannot be taken as independent ones, contrary to what has been 
done in the literature. By using these relations we see that 

= 0. (185) 

Note that this condition can be obtained by imposing just IR regularity conditions, and therefore 
any solution with a non-zero is singular in the IR. Imposing all those conditions, the 1/r 4 terms 
agree with those of [14]. However, the agreement is not complete: there is some discrepancy 
with the non-logarithmic term in <pf. 
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